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Abstract 

We investigate the current noise of a charge qubit coupled to a phonon bath in different parameter regimes. We find, using 
the theory of Full Counting Statistics of non-Markovian systems, that the current fluctuations are strongly influenced by 
memory effects generated from the interplay between quantum coherence and the dynamics of the phonon bath. 
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The recent experimental realizations of charge 
qubits in double quantum dots clearly illustrate the 
importance of controlling quantum coherence within 
a system [1]. Detecting quantum coherence via indi- 
rect and non-invasive measurements is thus becoming 
a crucial task. Currents fluctuations, in particular 
the current noise, have been suggested as a useful 
investigation tool [2,3]. In a recent experiment, fluctu- 
ations of the current through two coherently coupled 
quantum dots were measured [4], and the observed, 
unexpected super-poissonian noise was subsequently 
explained as a pure effect of quantum coherence [5] . In 
the present work, we consider the same system, how- 
ever, using an approach based on the theory of Full 
Counting Statistics for systems with Non-Markovian 
Dynamics (FCSNMD) [6,7]. The main advantage of 
this technique lies in the possibility to calculate all or- 
ders of the current cumulants using a unified approach 

^ Corresponding author. E-mail: braggio@fisica.unige.it 



that allows us to treat not only the Markovian, but 
also the non-Markovian dynamics of the charge trans- 
port process. We show that the observed, asymmetric 
super-poissonian behavior is determined by the inter- 
play between memory effects induced by the coherent 
dynamics and the coupling to the dissipative phonon 
bath. This indicates that it may be disputable to asso- 
ciate super-poissonian noise purely to coherent effects. 
We conclude that current fluctuations in general can 
be strongly modified by memory effects that may be of 
both intrinsic and/or extrinsic origin to the system [7]. 

We consider the model of transport through a charge 
qubit in a dissipative environment previously studied 
in Refs. [2,3]. The model consists of a double quantum 
dot operated in the strong Coulomb blockade regime 
and tuned to a degeneracy point, where electron trans- 
port proceeds only via the three charge states |0) (un- 
occupied dots), |+)/|— ) (left/right dot occupied). The 
Hamiltonian of the system reads 
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(1) 



where we denote the tunnel coupling between |+) and 
|— ) by Tc and the energy difference (detuning) of the 
two states by e, having introduced the standard spin- 
boson pseudo-spin operators CTz = |+)(+| — K )(— I 
and ffx = l+X^ 1 + Relaxation and dephasing 

are assumed to occur due to the coupling Vbo-z = 
gaz X^jC'^l + o.j)/'^ to a surrounding dissipative heat 
bath described as a reservoir of non-interacting bosons, 
Hb = X]j fit^JjOjCLj. We are thus dealing with the well- 
known spin-boson problem. Finally, we introduce a 
left/right (a — +, — ) lead described as non-interacting 



fermions, 



Ck,a coupled 



to the spin-boson system via the tunnel-Hamiltonian 
^T = Efe,„=±(iccIjOXa|+h.c.). 

To describe charge transport through the system we 
follow the approach developed by Gurvitz and Prager 
[8]. Assuming a constant tunneling density of states 
and a large bias across the double-dot, we can derive 
an equation of motion for the reduced density matrix 
p'"' , defined in the Hilbert space of the electronic sys- 
tem (|0), j-f), j— )) and the phonon bath, resolved with 
respect to the number of electrons n which have tun- 
neled from the left barrier [8] . To calculate the current 
cumulants it is convenient to introduce the counting 
field X Buud consider the Fourier transform /S*--^^ — 
E^ e'"-^/?'"' [9]. Current moments can be obtained 
from the dynamics of the sum of the ^-dependent oc- 
cupation probabilities contained in the vector P{x,t) 
with elements [P{x,t)]i = {*|Trs{p('^' i = 
0, + , — . Here, Tr_B denotes a partial trace over the 
bath degrees of freedom and it is easy to verify the 
normalization for x = 0. We proceed by assuming a 
local Born approximation for the diagonal parts of 
the density matrix: at any time the bath is assumed 
quickly to reach local equilibrium corresponding to 
the given charge state. Within this approximation 
we may solve the problem for small values of TcO] 
The rate equation for the counting probability vector 
is P{x,t) = dT'W{x,t — t)F{x,t) with the mem- 
ory kernel W(x, t) containing the rates for tunneling 
to and from the leads and phonon-assisted rates for 
tunneling between the two dots. In the following we 
apply some recently developed results for the general 
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theory of the FCSNMD to calculate the current and 
the noise [7]. The starting point is the non-Markovian 
equation of motion for P{x,t) which in Laplace space 
translates to [z - W(x, z)]P{x, z) = P(x, t = Ojwith 
P(x,t = 0) being the initial condition aX t = OLH For 
the problem at hand, the memory kernel reads 



W(x,^) = 
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where r± = 27r£i±|t±| are the lead tunneling 
rates for energy independent densities of states 



Q±. The phonon assisted hopping rates are Pi*'' — 
T^[g'^^\z+) + where z± = z + T-/2 ± ie 

and the bath correlation functions in Laplace 
space are defined as g^'^\z) = /(,°§te-^*C;(±'(t). 
In time domain these functions read = 
Trs{e"'^B^'V^'^'e'^B~''} with the equilibrium den- 



/TrB{e-'^"B } of the 



r{±) 



fit) —[iH^ 

sity matrix = e ^ b 

displaced phonon bath with Hamiltonians H'^g — 
Hb ± Vb and inverse temperature /3 = I/ZcbT. Assum- 
ing an Ohmic spectral density, J{lo) = 2g^(x>e~^^^'' , 
we obtain to order g^ in the bath coupling and to 
leading order in 1/ /^ojc , the expression 



z 



27r / \ liJc 



(3) 

where H(x) = ln(a;) - *(a;) - l/2a; with be- 
ing the digamma function and V^'^^ix) = [sin(a;) =p 
i cos(a;)] [7r/2—si(a::)=i=ici(a;)], where si(a;) = dtsm{t)/t 
and ci(a;) = — dtcos{t)/t. We note that even with- 
out the phonon bath, corresponding to g — 0, we still 
have a contribution of the form g^^\z) = 1/z oc- 
curring due to the coherent dynamics, and only if we 
take the z — > limit we recover the standard Fermi's 
golden rule rates for incoherent tunneling between the 
dots, r^t^'o = T^r-/[{r-/2f + The current / and 
the noise Fano factor F = S{0) /el for the system can 
be obtained using the theory of FCSNMD [6,7] [3 The 
theory shows that noise and higher order cumulants, 
but not the mean current, may contain signatures of 



to all orders in Tc [2] . 



^ For the zero-frequency fluctuations considered here, the 
initial conditions do not play a role. 

* The theory also allows us to calculate higher cumulants, 
which we, however, will not consider here. 
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Figure 1. Fano factor vs. level detuning e for different phys- 
ical parameter regimes and two temperatures T = 1.4 K 
(light grey) and T = 12 K (dark grey). The other param- 
eters correspond to those used in Fig. 2b of Ref. [5]. Co- 
herent regime (solid black line), g = sequential tunneling 
(black dashed line), phonon-assisted incoherent tunnelling 
(dashed grey lines), phonon-assisted tunneling (solid grey 
lines) and, for comparison, perturbative approach in the 
hybridized basis (dotted lines). 

non-Markovian dynamics, signaled by the presence of 
z-derivatives of the rates ri*'. The current is given as 

/ = r_r+r<+Vr? with r? = r_(r+ + ri+^) + r+r^ 

and Fp = (r[,~' -I- r[,+^). The Fano factor F reads 

i' — fjn H p4 , 

^ 2r-r^r^-> + v\vl + r^_[r^ + [v^+^f] 

J- t 

(4) 

where Fm denotes the Markovian contribution to the 
Fano factor (taking only the 2 = contribution) and 
where we have defined SzT^*-' s 9z[ri^']2=o. 

In Fig. 1 we show the Fano factor calculated in 
different parameter regimes. Without coupling to the 
phonons (5 = 0), we find in the limit z — > the sequen- 
tial tunneling regime described only by the Markovian 
terms [6]. The Fano factor is symmetric around e = 
and sub-poissonian (black dashed curve). Taking the 
full z-dependence of the rates, non-Markovian con- 
tributions are included, and we find the so-called 
coherent regime [8,10]. The Fano factor is again sym- 
metric, but now super-poissoniam (solid black curve) 
[5]. Including coupling to the phonon bath {g > 0), 
we are in the phonon-assisted regime, where electron 
tunneling between the dots is assisted by phonon 
absorption and/or emission. In the incoherent ap- 
proximation without non-Markovian contributions, 
we find sub-poissonian noise independently of the 



temperature (dark and light grey dashed lines). On 

the other hand, including non-Markovian terms we 
see the asymmetric super-poissonian behavior for low 
temperatures (dark grey solid lines) which disappears 
with increasing temperature (light grey solid lines) as 
also seen in the experiment [4]. We, moreover, find an 
optimal agreement with the results of the standard 
perturbation theory in the hybridized basis (dark and 
light grey dotted lines), allowing us to exclude possi- 
ble spurious effect of the perturbation in Tc. We see 
that the noise becomes super-poissonian duo to the 
non-Markovian corrections. By direct inspection of 
the expression for the Fano factor, we see that the 
condition for having super-poissonian noise is given by 
the inequality r+(9,r^+'(r_ + v'f'') - Fj^+^a^F^"') > 
rQ'''^(rp -h r_ -l- r+). This shows that a pure coherent 
contribution cannot easily be disentangled from other 
classical or quantum bath memory effects, since these 
all appear mixed in the memory kernel. Our analy- 
sis supports the suspicion, previously arisen, that an 
anomalous behavior of the current cumulants may 
be due to classical memory effects [7]. Such effects, 
independently of their origin, can destroy the signa- 
tures of pure quantum efi^ects, and should be carefully 
taken into account in any interpretation of current 
fluctuations. 
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